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ABSTRACT 

Aims. The aim of this work is to study the energy transport by means of MHD waves propagating in quiet Sun magnetic 
topology from layers below the surface to the corona. Upward propagating waves find obstacles, such as the equipar- 
tition layer with plasma /3 = 1, the transition region and null points, and get transmitted, converted, reflected and 
refracted. Understanding the mechanisms by which MHD waves can reach the corona can give us information about 
the solar atmosphere and the magnetic structures. 

Methods. We carry out two-dimensional numerical simulations of wave propagation in a magnetic field structure that 
consists of two vertical flux tubes with the same polarity separated by an arcade shaped magnetic field. This configu¬ 
ration contains a null point in the corona, that significantly modifies the behaviour of the waves as they pass near it. 
Results. We describe in detail the wave propagation through the atmosphere under different driving conditions. We also 
present the spatial distribution of the mean acoustic and magnetic energy fluxes for the cases in which these calculations 
are possible and the spatial distribution of the dominant frequencies in the whole domain. 

Conclusions. We conclude that the energy reaches the corona preferably along almost vertical magnetic fields, that is, 
inside the vertical flux tubes. This energy has an acoustic nature. Most of the magnetic energy keeps concentrated 
below the transition region due to the refraction of the magnetic waves and the continuous conversion of acoustic-like 
waves into fast magnetic waves in the equipartition layer located in the photosphere where plasma /3 = 1. However, 
part of the magnetic energy reaches the low corona when propagating in the region where the arcades are located, but 
waves are sent back downwards to the lower atmosphere at the null point surroundings. This phenomenon, together 
with the reflection and refraction of waves in the TR and the lower turning point, act as a re-feeding of the atmosphere, 
that keeps oscillating during all the simulation time even if a driver with a single pulse was used as initial perturbation. 
In the frequency distribution, we find that high frequency waves can reach the corona outside the vertical flux tubes. 
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1. Introduction 


Quiet sun magnetic structures are much less organized 
and have less intense magnetic field compared to sunspots 
(Sanchez Almeida & Martinez Gonzalez 2011). At the same 
time, quiet Sun covers at least about 90% of the solar sur- 
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magnetic structures to be important contributors to the 
solar magnetism. These weak structures are very dynamic, 
con stantly perturbed by convection and propagating waves 
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tion at the base of the photosphere (Goldreich fc Keeley 

1977 Balmforth 1992 Nordlund &: Stein[ ~ 2QQip . TEe" mag¬ 

netic held concentrations embedded in the intergranular 
lanes suffer foot point motions and this d rives waves in these 


structures throug h the solar atmosphere (Hasan et al. 2000 
Kato et al. 2011). Depending on the foot point motion, dif- 


ferent waves are going to be excited. There is a large num¬ 
ber of works on numerical simulations that reproduce these 
different scenaria: horizontal motions producing transverse 


kink waves; pressure fluctuations 
waves; twisting motions generating 
observationally-driven or impulsive 
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some exceptions (]Choudhuri et al.| 1993 De Pontieu et al. 

2005; Fedun et al.| 2011a|), these works model the wave 


propagation in the photospheric and chromospheric layers. 


The dominant oscillations observed in the quiet Sun 
photosp here are acoust i c-gravity waves wit h a period of 5- 


minute (Leighton et al. 1962 Ulrich 1970). As these waves 


propagate upwards the frequency distribution changes. 
Essentially 3-minute period waves are observed in the 
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lation s appear in the corona (De Moortel & Nakariakov 
2012). The distribution of the dominant frequency of oscil- 


flux reaching the corona, as well as the frequency distribu¬ 
tion of waves with height and horizontal distance. 


2. Numerical Method 

We solve the two-dimensional ideal magnetohydrodynamic 
equations of conservation of mass, momentum, energy, and 
the induction equation for the magnetic field: 


lations along the atmosphere appears to be different inside 
the quiet network cells and network borders and facular 
regions, containing stronger magnetic fields. There are ob¬ 
servations of the 5-min waves channelled along the inclined 
magnetic field lines to the chromosphere in facular regions 
(via so-called “ramp-effect” or “magnetoacoustic-portals”), 
which have been explained by the reduction of the acoustic 
cut-off frequency in inclined magnetic fields (Michalitsanos[ 
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period oscillations vertically upward (see Roberts 
The latter mechanism, how ever, was not found to oe doni- 
inant in the simulations by Heggland et al. (2011), with a 
complex treatment of radiative transfer. 

As the waves propagate upwards they find obstacles, 
such as the equipartition layer where plasma /3 = 1 , the 
layer where the local acoustic cut-off frequency is equal to 
the wave frequency, the steep temperature gradient at the 
transition region, etc. In these regions waves suffer transfor¬ 
mations, refraction, and reflection. Since the atmospheric 
properties change as we move from one layer to another, 
and even from one column to another, these effects can 
appear at different horizontal and vertical positions. For 
instance while in the photosphere the gas pressure domi¬ 
nates, the magnetic pressure becomes dominant from the 
middle chromosphere upwards. 

Particularly interesting is the question of how, and 
in which amount, the energy of waves excited in sub- 
photospheric layers can reach the chromosphere and corona, 
and whether the wave modes reaching there can be dissi¬ 
pated to convert their energy into heat, i.e. how the so¬ 
lar atmospheric layers are magnetically and energetically 
connected by means of waves. Only few numerical studies 
perform simulations of waves in the whole domain from 
the photosp here to the chromo sphere, transition region 
and corona. Fedun et al. (2011a), in their 2D simulations 
of high-frequency wave propagation, find that magneto¬ 
acoustic waves in a wide range of frequencies can effec¬ 
tively leak energy into the corona. These authors consid¬ 
ered an isolated magnetic fiux tube that becomes vertical 
in the chromosphere. In the present work we study how 
the energy reaches the chromosphere and corona by excit¬ 
ing waves below the photosphere with different perturba¬ 
tions. We choose a complex magnetic field that consists 
of vertical flux tubes of the same polarity separated by an 
arcade-shaped magnetic field, resembling a network and in¬ 
ternetwork regions, with wave periodicities in the 3-5 min 
regime. We discuss in detail the behaviour of different wave 
modes in such complex situation. The length of our simu¬ 
lated time series allows us to reach the stationary regime in 
two of the three considered cases and quantify the energy 


where all the notations are standard. We use an ideal 
equation of state and 7 = 5/3. The terms with subscript 
“diff” are the artificial hyper-diffus ive terms req u ired fo r 
the code’s numerical stability, see Felipe et al. l (|2010a). 


After removing the equilibrium condition (see below), we 
solve these non-linear equations by means of the code MAN¬ 
CHA w hich is described in deta il in IKhomenko fc Collados] 
(2006) and Felipe et al. (2010a). The numerical code solves 
the equations for perturbations. Through the present study 
we keep the amplitude of the perturbations small enough 
to concentrate on the linear regime. The non-linear effects 
will be reported in a forthcoming study. 

2.1. MHS model 

For simplicity, we have chosen a potential magnetic field 
structure. Such magnetic field does not interact with the 
equilibrium atmosphere. This allows us to compute the hy¬ 
drostatic equilibrium model separately from the magneto¬ 
static equilibrium. 


2.1.1. Hydrostatic model 

To get the hydrostatic equilibrium we start by defining the 
temperature structure as a function of the vertical coor¬ 
dinate 2 : = {—5,10} Mm. For that, we join three models. 
The sub-photospheric la yers are described by the convec - 


tively stabilized model byfParchevsky & Kosovichev (2007). 
The main reason to choose this model is because no con¬ 
vective instability is developed in response to the pertur- 



mal corona with a temperature of one million Kelvin. The 
complete temperature distribution is shown in Figure 
(green line). Once the temperature stratification is estab¬ 
lished we calculate the scale height i7, and then integrate 
the equation for pressure and recover the density distribu¬ 
tions, as follows: 


H = 


RgasTp 
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dpojz) ^ 
dz H 



(6) 

( 7 ) 


using the ideal gas equation. With this, we ensure that 
the atmosphere is in hydrostatic equilibrium. The result¬ 
ing distributions of pressure and density are given in Fig. 

A plane-parallel atmosphere is built by replicating the 
obtained vertical stratifications. In the above equations we 
denote the equilibrium variables by a subscript “ 0 ”; Rgas 
is the gas constant and /i is the mean atomic weight. We 
set fi=l at the lower layers of the atmosphere where we 
can consider that the gas is made of neutral hydrogen, and 
/i=l /2 at higher layers where the plasma is fully ionized. 


2.1.2. Magnetostatic model 

The magnetostatic model is given by an arcade shaped mag¬ 
netic field plus a constant component: 



0 2 4 6 8 10 12 14 

HORIZONTAL [Mm] 


Bqx = Booex.p{—kz)sm{kx) ( 8 ) 

Boz = Bu ^ Booex.p{-kz)cos{kx) (9) 

where Bqo and are the magnetic field strengths in the 
photosphere and corona, respectively, k defines the spatial 
(horizontal and vertical) scales and x and ^ are the hori¬ 
zontal and vertical directions, respectively. We choose Bqo 
= 100 G and B^i=10 G, so we are working with a weak 
magnetic field region. Equations defines a magnetic field 
configuration composed by two vertical flux tubes separated 
by an arcade-shaped magnetic field, as drawn in Figure 
Ghoosing higher values for the photospheric and coronal 
magnetic field has the disadvantage of significantly reduc¬ 
ing the computational time steps due to the high values of 
the Alfven speed in the computational domain. Higher val¬ 
ues of Bqo and Bu would lower the eventual position of the 
l3 = 1 contour, but would not change the overall picture of 
the wave propagation and conversion. Table gives sound 
and Alfven speeds for the different layers in the equilibrium 
atmosphere. 

Our field configuration includes a null point in the 
corona. In this point the magnetic field strength is mathe¬ 
matically zero and the plasma (3 goes to infinity. This null 
point changes drastically the behaviour of the waves as they 
propagate near it. As the magnetic field strength is zero, 
the Alfven speed is null, and therefore, only the acoustic 
waves can physically pass through it. 

2 . 2 . Simulations setup 

We analyse three simulation runs, performed with the same 
model for the equilibrium atmosphere, but with a different 
driving of the waves. The cases include: (1) vertical driving, 
when the analytical solution for an acoustic wave with a pe¬ 
riod of 200 sec (fast mode) is applied at the lower boundary 
of the simulation domain; ( 2 ) horizontal driving, when the 
analytical solution for the slow magnetic wave with a period 
of 300 sec is applied at the lower boundary; (3) instanta¬ 
neous pressure pulse below the photosphere. Table gives 
a summary for these runs, including the section and the 
figures where the results are presented. We choose the peri¬ 
odicities of 200 and 300 sec to study the behaviour of waves 
with frequencies below and above the cut-off frequency, typ¬ 
ical for the solar conditions. The cut-off frequency layer 


Fig. 2. Magnetic field strength distribution from below the 
photosphere (—5 Mm) to the corona (10 Mm). The vertical 
axis gives the height in Mm and the horizontal axis gives 
the horizontal size in Mm. Height z=0 represents the base of 
the photosphere. The red lines are the magnetic field lines 
and the yellow dashed lines are the [3 =l contours (one on 
the photosphere and the other one around the null point). 


for 300 sec period waves is located below the photosphere, 
those for 200 sec waves is at the temperature minimum. 


2.3. Boundary conditions 

To avoid spurious reflections of waves at the upper bound¬ 
aries of the domain, w e add a Perfectly Ma tched Layers 


119941) of 20 grid points, see Felipe et al. 
ils. The PML is specially designed to ab- 


(PML [B erenger 
j 2010 bp for detai 
sorb the perturbations that reach the borders of the do- 


tions of MHD waves (see 

Parchevsky & Kosovichev, 

2008 

Khomenko et al.[|2008t|Fe 

iipe et al.||2010bt|Hanasoge et al. 

2010). In the simulations with a periodic driver, the lower 
boundary is replaced by the analytical solution for a given 


type of wave. In the simulations with an instantaneous 
pulse, the lower boundary contains the PML layer as well. 
In the horizontal direction, boundary conditions for peri¬ 
odic variations are used. 


3. Driving waves with linear perturbations 

This section presents the results of the different simulation 
runs. Despite the code can handle non-linearities, the am¬ 
plitude of the perturbation is sufficiently small and the wave 
propagation is kept in the linear regime in all the cases, so 
no shocks are formed. This is done on purpose because our 
intention is to study the wave propagation, conversion, re¬ 
fraction and reflection in the linear regime. To better under¬ 
stand the results of the simulations, we briefly summarize 
the nomenclature, and some basics of the mode conversion 
theory. 
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Fig. 1. Temperature (green line), density (blue line) and pressure (red line) distributions. The left vertical axis gives us 
the values for temperature and pressure, and the vertical axis on the right is the density scale. 


Table 1. Alfven and sound speeds in the different layers of the equilibrium atmosphere. 


Photosphere 

Chromosphere 

Corona 

Cso 8 km/s 

Vao 2 km/s 

10 km/s 

100 km/s 

100 km/s 
1000 km/s 


Table 2. Summary of the simulation runs 


Section 

Driving 

Amplitude 

Period 

dx / dz /Duration 

Figures 

3.1 

vertical, harmonic 


200 s 

75 km / 50 km / 3000 s 

3-6 

3.2 

horizontal, harmonic 

0.2 m 

300 s 

75 km / 50 km / 3000 s 

7-10 

3.3 

pressure pulse 

lO-^'po 

- 

75 km / 50 km / 3000 s 

11-13 


3.1. Wave transformation 

Dashed-dotted curves in Figure show the f3 = 1 contours 
in our equilibrium model, one in the photosphere and an¬ 
other one around the null point. The plasma (3 gives us 
the information of whether the atmosphere is magnetically 
dominated or not. Strictly speaking, the plasma /3 is defined 
as /3 = Pgas/^mag 5 but from the point of view of the wave 
propagation, it is more convenient to define it as a squared 
ratio of the characteristic speeds, i.e. /3 = (?slv\. The dif¬ 
ference between both definitions is a factor 7 / 2 , which in 
our case is given by 5/6. This value is close to unity and, 
therefore, almost imperceptible in the figure or in wave ef¬ 
fects. We will use the latter definition through the paper. In 
the layer with /3 ~ 1 (equipartition layer), transformation 
of different wave modes can occur. Since our simulations 
are 2D, the only possible wave modes are fast and slow 
magneto-acoustic-gravity modes, so we exclude the conver¬ 
sion to Alfven mode from our current analysis. 

The propagation of the fast and slow magneto-acoustic- 
gravity modes is in general rather complex, except in cases 
with very simple magnetic topologies. Under a mathemat¬ 


ical point of view, the analysis of the behaviour of waves 
comes in a natural way if done according to the fast or slow 
nature of their propagation speed, regardless of their acous¬ 
tic or magnetic properties. However, for the interpretation 
of observations, it is in general more convenient to take into 
account the physical nature of the main restoring force. It 
is unfortunate that the terms wave transmission and con¬ 
version have different meanings in these two descriptions. 
Since this paper focusses on the interpretation of numer¬ 
ical simulations, we have preferred to follow the physical 
convention. This way, we will consider that a wave is trans¬ 
mitted when the nature of its main restoring force does 
not change when crossing a given layer, keeping its mainly 
acoustic/magnetic nature (independently of whether in the 
mathematical description it may change from fast to slow 
or viceversa). Similarly, we will use the term wave conver¬ 
sion when the nature of the main restoring force changes 
from acoustic/magnetic to magnetic/acoustic (keeping the 
fast or slow character or the wave). 

Depending on the magnetic field inclinatio n, the conver¬ 
sion and transmission are usually only partial (Zhugzhda fc 


Dzhalilovl [T9821 [Cjdf^ [20051 [20061 ISchnnker fc Cally[ |2006| 
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Khomenko & Collados 2006 Khomenko & Cally 2012 ), 
hence part of the wave keeps its acoustic/magnetic nature 
and part of it changes. The transmission efficiency depends 
on the attack angle, i.e., the angle between the m agnetic 


field lines and the wave vector in the layer of /3 = 1 (Cally 
200 ^: 


T = exp 


Kn sin^ 


{d/ds){cl/vl) 


( 10 ) 


where a is the attack angle and s is the distance along the 
direction of k and K = \k\. Cg/va is the ratio of the speed 
of sound and the Alfven speed. The transmission is com¬ 
plete when the magnetic field and the direction of the wave 
propagation are aligned, a = 0. In this case, the conversion 
coefficient, 1 — T, is zero. When a 7 ^ 0, then the transmis¬ 
sion decreases with increasing k, that is, with increasing 
wave frequency. 


3.2. Vertical periodic driver of 200 seconds 


In this run we drive waves by applying a vertical periodic 
perturbation in velocity and thermodynamic variables at 
the bottom of the domain and let the perturbation prop¬ 
agate a total duration of 3000 seconds. This duration is 
long enough for the waves to fill the full domain and reach 
a stationary motion. The direction of propagation of the 
perturbation in this case is nearly parallel to the magnetic 
field lines, so acoustic like fast mode waves are mostly gen¬ 
erated. We use an analytical solution for an acoustic-gravity 
wave i n an isothermal atmosphere, see |Mihalas fc: Mihalas 
(1984). The perturbations are given by: 


Svz 

6p 

Po 

Sp 

Po 


Vo exp + k^izj sin(wi - k^rz) (11) 

Vo\P\ exp + k^iZ^ sin(wi - k^rZ + 4>p) ( 12 ) 

Vb|-R| exp + kzi^ sin(wt - k^rZ + 4>r) (13) 


where Vb = 10 “^ m s“^ is the amplitude of the velocity, see 
Table and H is the pressure scale height. Perturbations 
in the horizontal velocity and the magnetic field are null. 
Vzr and \izi are the real and complex parts of the vertical 
wave number,respectively. Given the wave frequency cj, this 
is found from the dispersion relation: 


— kzr “ 1 “ ikzi — ^ —- ( 14 ) 

where = ^gj^Cg is the isothermal acoustic cut-off fre¬ 
quency. 

The relative amplitudes and phase shifts between the 
perturbations are given by 


1^1 - j ATAvyA 
1 ^ 1 = 


4>p 


arctan 

arctan 



1 (7-2) \ 

2Hkzr 7 J 


2Hkzr 


( 15 ) 

(16) 

(17) 

(18) 


The perturbation is imposed in the interval z = 
{—5, —4.65} Mm, i.e. occupies 7 grid points. 

As is frequently done in MHD wave studies, to distin¬ 
guish between the modes in the simulations, we project the 
velocities into the directions parallel and perpendicular to 
the local magnetic field, defined as: 

^long = CO^OOz (19) 

etrans = - ^ ^vnOcz ( 20 ) 

where 0 is the inclination of the magnetic field. This way, 
in the ^ 1 atmosphere, the acoustic slow mode is visi¬ 

ble in the longitudinal velocity (propagation parallel to the 
magnetic field); and the fast magnetic mode is essentially 
transverse to the magnetic field. Note that the projections 
do not help to distinguish between the modes in the ^ 1 

region. 


3.2.1. Wave propagation 

Figure]^ shows the longitudinal (left panels) and transverse 
(right panels) velocity projections, scaled with the factors 
^^poCs and respectively. Such scaling gives an ap¬ 

proximation for the energies contained in the different wave 
components. Some important layers are also indicated in 
each panel of this figure: ^=1 contours (dashed black lines), 
the transition region (red line) and the location of the cut¬ 
off layer with a period of 200 seconds (black solid line), the 
same period as that of the driver. 

The magneto-acoustic fast and slow waves start propa¬ 
gating upwards through the layers below the photosphere, 
where plasma ^ 1 (see the rows t = 210 s and t = 360 

s). As the driver is parallel to the magnetic field at the foot 
points, the amplitude of the longitudinal velocity is larger 
than the transverse one. 

As the waves reach the equipartition layer at t = 450 
s, fast acoustic-like waves are partially converted into fast 
magnetic-like waves and partially transmitted into longitu¬ 
dinal slow acoustic waves. These can be seen comparing the 
second (360 s) and third (450 s) rows of Fig. in the sec¬ 
ond row the transverse component crosses the equipartition 
layer with a very small amplitude, but once the longitudi¬ 
nal component reaches the layer where /3 = I at 450 s, the 
amplitude of the transverse component grows considerably 
while the amplitude of the longitudinal waves decreases. 

These waves continue propagating upwards until they 
reach the transition region (before 600 s). At this point, 
part of the wave energy is reflected and part is transmitted 
into the corona. The transverse velocity component suffers 
significant reflection from the transition region (right pan¬ 
els from 750 s to the end). The energy of the transverse 
component is mostly concentrated between the equiparti¬ 
tion layer and the transition region. 

Some of the energy of the longitudinal waves is trans¬ 
mitted to the corona through the transition region. The 
behaviour of these waves, once they pass the transition 
region, depends on whether they propagate near the null 
point or far from it (see left panels of Fig. 1^ from t=750 
s to the end). The waves propagating near the lateral bor¬ 
ders of the domain, inside the almost vertical flux tubes, 
are longitudinal and continue propagating upwards along 
the field lines up to the upper corona. Waves propagating 
outside the vertical flux tubes reach the null point and this 
changes drastically their behaviour. These waves suffer par¬ 
tial conversion, from magnetic to acoustic and vice versa. 
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-400 -200 0 200 400 -200 -100 0 100 200 


V erg/(cm^ s) V erg/(cm^ s) 

Fig. 3. Longitudinal (left) and transverse (right) velocities 
multiplied by and respectively, for the sim¬ 

ulation run with a vertical periodic driver. The black solid 
line is the location of the cut-off frequency for the 200 
second period waves, the dashed black lines are the /3=1 
contours and the red line shows the transition region. The 
magnetic field lines are inclined black lines. Positive veloc¬ 
ity (red color) is an upflow. Note that the horizontal vari¬ 
ations in the longitudinal and transverse velocities in the 
bottom layers where the driver is imposed, despite the ve¬ 
locity perturbation is constant in the horizontal direction 
(see Eq. pTj ). These changes are caused by the projection 
of the velocity on the directions parallel and perpendicular 
to the magnetic field, which has a varying inclination. The 
movie of the wave propagation is available in the online 
version of the paper. 
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Fig. 4. Time-height diagram of y^po^'^iong (left) and 
x/Po'^a'^trans (right) for three selected horizontal positions in 
the vertical periodic driver simulations. Upper row: x = 150 
km, where the magnetic field is nearly vertical. Middle row: 
X = 3750 km, close to the center of the left arcade with 
nearly horizontal magnetic field. Bottom row: x = 7500 
km close to the null point location. The color scale is the 
same as in Figure We notice how the simulation reaches 
a stationary stage after about 1000 s. 


and transmission, from fast to slow and vice versa, at the 
equipartition layer around the null point. The Alfven speed 
is zero at the null, and therefore, only acoustic fast waves 
(inside the circle marked by the coronal /d = 1 contour) 
are able to cross it. The movie of the simulations available 
in the online version of the paper) makes apparent how 
the null point “absorbs” the waves and then “ejects” them 
away in all directions. This effect is mainly caused by the 
refraction of the fast magnetic waves due to the large gradi¬ 
ent of the Alfven speed. Part of these fast waves propagate 
downward again through the transition region and equipar¬ 
tition layer to the lower atmosphere. The null point acts as 
a re-feeding of the atmosphere. 

Few minutes after the waves reach the upper bound¬ 
ary, the simulation enters into the stationary stage. This 
property will be used for the calculation of the mean wave 
energy fluxes, done in the section below. 

Time-height diagrams allow to appreciate more clearly 
the behaviour of waves described above, see FigureThere 
we plot the quantities y^po^Uong (left) and ^po'^a'^trans 
(right) at three different horizontal locations in the domain. 
The cuts are located at x = 150 km inside the flux tube. 
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where the magnetic field is nearly vertical, at x = 3750 
km in the center of the left arcade and where the magnetic 
field is nearly horizontal, and at x = 7500 km, in the middle 
of the domain where the null point is located. The inclina¬ 
tion of ridges on this diagram reflects the wave propagation 
speed (more vertical means faster). Fig. 3 shows how the 
picture of wave interference changes at dmerent magnetic 
field inclinations. 

The left panels demonstrate that waves propagating 
parallel to the nearly vertical magnetic field at x = 150 km 
and at X = 7500 km travel almost undisturbed up to the 
transition region, where reflection is produced and an inter¬ 
ference pattern is formed between the upward and down¬ 
ward propagating waves. As the attack angle is small at the 
/3 = 1 contour at these locations, the transmission from fast 
acoustic to slow longitudinal acoustic waves is almost total. 
Once the waves at x = 150 km reach the transition region 
at about 600 s of the simulation time, they continue prop¬ 
agating with significant amplitude upwards (top left panel 
of Fig. [3. Unlike that, waves propagating at x = 7500 km 
are affected by the null point. One may appreciate at the 
bottom left panel, how the waves propagating downward, 
due to the refraction suffered near the null point, meet the 
waves going upwards and their amplitudes cancel out. 

The behaviour of waves propagating at x = 3750 km 
(middle panels of Fig.|^, where the magnetic field is nearly 
horizontal below the transition region, is different from 
above. One can observe several small-scale wave compo¬ 
nents, with their energy located essentially between the 
/3 = 1 contour and the transition region (middle left panel). 
These small-scale waves are a result of the mode conver¬ 
sion into slow modes propagating along the inclined field 
(their propagation speed is apparently small because of the 
longer distance they have to cross in the inclined fields). 
Fast magnetic modes have much larger propagation speeds 
and can be seen at the middle right panel as nearly verti¬ 
cal lines with a maximum power concentrated again in the 
chromosphere below the transition region. The reflection of 
longitudinal slow waves from the null point is also apparent 
at the middle left panel from the opposite inclination of the 
ridges above the transition region and below the null point. 


3.2.2. Energy fluxes 

We calculate the mean acoustic and magnetic energy fluxes 
as follows (Bray & Loughhead, 1974): 

i^ac = {SpSv), (21) 

-F’mag = {SB X {5v X Bo))/lJo . (22) 

In order to obtain the meaningful averages of the flux one 
needs to include at least several periods of the stationary 
regime of the simulations. For the simulations presented 
above, we considered that the stationary regime is main¬ 
tained from 1000 to 3000 s of the simulation time. 

Figure shows the spatial distribution in the domain 
of the magnetic and acoustic, vertical and horizontal mean 
energy fluxes. It can be seen that almost all the magnetic 
energy (lower panels) is concentrated between the photo¬ 
sphere {z = 0 km) and transition region (z = 2100 km), 
while the upward acoustic energy flux is significant in the 
corona inside the vertical flux tubes. The only significant 
upward energy flux in the corona is due to upward prop¬ 
agating longitudinal slow acoustic waves. No wave energy 
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Fig. 5. Mean energy fluxes for the simulation with the verti¬ 
cal periodic driver. Top: horizontal and vertical mean acous¬ 
tic energy fluxes. Bottom: horizontal and vertical magnetic 
mean energy fluxes. The lower black line is the acoustic cut¬ 
off frequency layer for 200 second period, the dashed lines 
are the /3 = 1 contour, the red line is the transition region. 
The magnetic field lines are inclined black lines. The units 
of the bar are 10^ [erg/s/cm^]. 


flux is present around the null point. This is happening 
because waves propagate in all directions around the null 
point and cancel out when waves with opposite directions 
meet. 

In Figure three different cuts of the mean vertical en¬ 
ergy fluxes are shown in more detail. As before, we made 
cuts at three horizontal locations: inside the flux tubes at 
X = 150 km; at the middle of the arcades where the mag¬ 
netic field is horizontal at x = 3750 km and at x = 7500 km 
where the null point is located. The bottom panel of Fig.|^ 
clearly shows that all the magnetic energy is concentrated 
between the photosphere {z = 0 km) and the transition 
region, as a consequence of the mode conversion in the in¬ 
clined magnetic field of the arcades (green line). The sign of 
the energy flux is negative meaning that waves essentially 
propagate down, which is explained by the influence of the 
null point. The magnetic energy flux in our experiment is 
unable to cross the transition region. 

The upper panel of Fig. shows that the acoustic en¬ 
ergy in the photosphere propagates upwards in the nearly 
vertical magnetic concentrations (blue and red lines), and 
downwards at the locations of the arcade. This acoustic en¬ 
ergy flux is significantly larger in the photosphere and in 
the coronal part, compared to the magnetic flux (the latter 
is almost zero in these regions). The largest upward acous¬ 
tic energy flux reach the corona in the nearly vertical flux 
tubes at the sides of the domain and makes less than 10- 
20 % of the original energy flux of acoustic waves below the 
photosphere. 

3.3. Horizontal periodic driver of 300 seconds 

In this run, a horizontal periodic motion with a period of 
300 seconds is applied at the lower boundary of the simu- 
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Vertical Mean Acoustic Flux 




Height [Mm] 


Fig. 6. Vertical cut of the vertical component of the acous¬ 
tic (top) and magnetic (bottom) mean fluxes for three dif¬ 
ferent horizontal positions, as indicated in each panel. The 
vertical black lines are the locations of the cut-off frequency 
for 200 seconds period waves, /3 = 1 contours (dashed), 
transition region and the null point (from left to right). 


lation domain at heights z = {—5, —4.65} Mm, same as for 
the vertical driver. 

= Vq sin(cc;t) (23) 

^ = ’^Vosiniut) (24) 

no uj 

where the horizontal wavelength is infinite, so kx = 0, 
and kz = uj/va-The perturbation in the rest of the vari¬ 
ables is null. The value of the magnetic field strength 
Bo = ^/W. + is taken at the boundary. The horizon¬ 
tal velocity amplitude is chosen to be of Vo = 0.2 m s“^, 
see Tab. Such driving generates essentially no acoustic 
waves in vertical magnetic fields where the motion of the 
driver is completely transverse. This can be seen in some 


points of the edges of the domain and in the middle, be¬ 
tween the two arcades. However there appear acoustic-like 
waves in almost the rest of th e domain, where the mag¬ 


netic field is more incline d (see [Hasan et al. 
& van Ballegooijen 2008| ) 


2005 Hasan 


3.3.1. Wave propagation 

Figure shows the propagation of waves generated by the 
horizontal periodic driver. The driver excites magnetic slow 
waves seen as a small wavelength pattern at the bottom of 
the domain at all times. The pattern due to these waves 
disappears after some height around —3 Mm since these 
small-wavelength perturbations are affected by the numer¬ 
ical diffusion effects and are not fully resolved. At the loca¬ 
tions with inclined field, fast acoustic waves are also excited. 
These can be appreciated as a large-scale disturbance prop¬ 
agating from the bottom boundary upwards at the two up¬ 
per panes and at later times as well. These fast waves reach 
the equipartition layer at t = 300 s, and they get partially 
converted and transmitted. The perturbation seen in lon¬ 
gitudinal velocity (left panels) remains almost completely 
below the equipartition layer at all posterior time moments. 
The only exception are the vertical flux tubes at the domain 
sides where again, similar to the case of the vertical driving, 
some perturbation propagates upwards through the transi¬ 
tion region to the corona. 

Perturbations seen in the transverse velocity (right pan¬ 
els) above the equipartition layer and below the transition 
region at later times (t > 450 s) correspond to fast mag¬ 
netic waves. Neither these fast magnetic waves, nor slow 
acoustic waves (small-scale perturbations between the layer 
where /3 = 1 and the transition region at the left panels), 
pass through the transition region with large amplitudes, 
and their energy remains essentially below it. The ampli¬ 
tudes of waves reaching the corona are smaller than in the 
case of the vertical driver. 

Unlike the longitudinal slow waves in vertical flux 
tubes, fast magnetic transverse waves reach the corona in 
the whole domain. When passing through the null point, 
they are partially converted into acoustic-like waves in the 
equipartition layer around it. The simulation reaches a sta¬ 
tionary state after about t = 1200 seconds of the simula¬ 
tion. 

Figure [S] shows the time-height diagrams of 
(left) and -^/po'^a'^trans (right) at three different horizontal 
positions of the domain. The top and bottom rows of this 
figure (horizontal positions x = 150 km and x = 7500 km) 
show waves propagating in the nearly vertical magnetic 
field. One can appreciate how the behaviour of these waves 
changes at the equipartition layer (dashed-dotted line) and 
at the transition region. The influence of the cut-off layer is 
not so clear in these diagrams. The transverse velocity prop¬ 
agation at these two positions of the domain is very similar. 
The mode conversion results in an increase of the transverse 
wave component as the waves cross the equipartition layer 
in the downward direction. The ridges on the diagram show 
the opposite inclination indicating downward propagation 
produced by the reflection at the transition region. In the 
corona, the ridges are vertical, indicating that waves prop¬ 
agate with Alfven speed, so their nature is magnetic. As 
these fast waves cross the null point, they are converted at 
the equipartition layer around it (bottom right panel of Fig. 
[^. The propagation speed of the waves is altered, which re- 
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suits in a shift of the position of the ridges at heights around 
the null point. 

The longitudinal velocity behaviour at x = 150 km and 
at X = 7500 km is quite similar, except of the amplitude 
of the waves. Some weak longitudinal slow waves appear 
above the transition region at x = 150 km, generated by 
horizontally propagating perturbations in the corona. 

The middle panel of Fig. ^ shows the wave propagation 
at the location of the inclined field of the arcade. There, the 
transverse waves propagate undisturbed up to the transi¬ 
tion region where they are reflected. Conversion from fast 
acoustic into fast magnetic waves occur at the equiparti- 
tion level as can be seen from the smooth transition in the 
inclination of the ridges at that location. Some weak acous¬ 
tic slow waves are also generated there (middle left panel). 
These waves propagate up to the transition region and are 
reflected there and form an interference pattern around 2 
Mm height. 


3.3.2. Energy fluxes 

Same as for the simulations with a vertical driver we calcu¬ 
lated the mean acoustic and magnetic fluxes, see Eqs. 
averaging over the stationary state of the simulations, from 
t = 1200 s to t = 3000 s (6 periods of 300 s waves). 

In the case of the horizontal driver the energy fluxes are 
very small, see Figure This figure shows the spatial dis¬ 
tribution of the mean energy fluxes. The upward acoustic 
energy flux is essentially concentrated below the arcades 
at heights below the transition region (upper right panel) 
and does not reach the corona. Curiously, the average up¬ 
ward acoustic flux is negative inside the vertical tubes on 
the sides, opposite to the case of the vertical driving. This 
might be related to the fact that in the simulations with 
the horizontal driver these waves are excited by horizon¬ 
tally propagating disturbances in the corona after the waves 
are refracted downward around the null point. The vertical 
magnetic energy flux (bottom right panel) is mostly posi¬ 
tive and is much higher above the photosphere than below. 

Figure shows the vertical magnetic and acoustic 
fluxes at three different horizontal locations. This figure 
reveals that the largest acoustic energy flux in this simu¬ 
lation run is concentrated below the arcades (green line at 
the upper panel). The magnetic flux is very small except 
at heights around the equipartition region. The values of 
the flux are much lower than in the vertical driving case, 
despite the its larger driving amplitude. 


3.4. Instantaneous pressure pulse 

In this run we exerted a localized instantaneous force at 
the base of the arcade, in the middle of the domain, that 
consists of a pressure pulse of a gaussian shape, as follows: 


Spi 

Po 

Spi 

po 


( 


exp 
= A exp — ^ 


{X - Xq)^ 

{x - Xo)‘^ 
2 ^^ 


+ 

+ 


2 ( 7 ^ J\ 

2a2 J\ 


(25) 

(26) 


where A=10“^ is the relative amplitude given to the pulse, 
[xq^zq] = [7500,-3500] km are the coordinates at which 
the pulse is located and, = 1500 km and az = 1000 km 
are the widths of the two-dimensional Gaussian profile, see 


Tab. Such a pulse generates a superposition of waves in 
a broad range of frequencies. 

Unlike the simulations with the vertical and horizontal 
driving, we locate a PML boundary condition at both, top 
and bottom, boundaries of our simulation domain. 


3.4.1. Wave propagation 


The time evolution of the scaled longitudinal and trans¬ 
verse velocity components of the pulse simulation are given 
in Figure The pulse starts expanding in all the direc¬ 
tions (two upper panels). Like in the other runs, the waves 
get partially converted and transmitted when reaching the 
equipartition layer around t = 300 s. The acoustic fast com¬ 
ponent of the pulse is partially converted and transmitted 
into magnetic fast and acoustic slow waves, respectively. 

The acoustic slow waves continue propagating upward 
and reach the transition region before t = 450 s. Part of 
these waves is transmitted into the corona and they reach 
the null point, where they become trapped and are sent 
again in all directions, except those waves inside the vertical 
flux tubes that continue propaga ting to upper layers along 
the field lines (left panels of Fig. [Ty . 

Acoustic waves reflected from the transition region 
reach equipartition layer again on their downward prop¬ 
agation, and they suffer another partial transmission. Fast 
acoustic waves produced after this secondary transmission 
propagate down back to the lower boundary of the domain, 
where they get refracted due to the gradients of the acous¬ 
tic speed and are finally reflected up around —4 Mm at 
the lower turning point. Therefore, these waves start prop¬ 
agating upwards again (see left and right panels of Fig. 
at t = 750 s at the lower part of the domain) suffering 
the same phenomena as before. Thus, despite our driver is 
not periodic, the atmosphere keeps oscillating reaching a 
nearly stationary regime without a significant decrease of 
the wave power during all 3000 s of the simulation time 
(only waves reaching the upper and lower PML boundaries 
are removed from the system). The null point together with 
the transition region and the lower turning point below the 
photosphere act as a re-feeding of the atmosphere, and the 
waves are trapped between these layers. 

The variations of the transverse velocity (right panels of 
Fig. 11) reveal that fast magnetic waves are produced after 


the mode conversion around 300 s, and then reach the tran¬ 
sition region and are partially transmitted to the corona 
before t = 450 s. These fast magnetic waves get refracted 
and reflected due to the gradients of the Alfven speed, and 
cross again the equipartition layer on their way back. On 
their way back, the fast magnetic waves are partially trans¬ 
mitted into slow magnetic waves with small wavelengths, 
that propagate downwards and their energy gets lost due 
to the numerical diffusion. At the moments posterior to 
t = 750 s, no perturbations in the transverse velocity can 
be appreciated. 

Same as for the other cases. Figure shows the time- 
height diagrams of the scaled longitudinal and transverse 
velocities at three representative horizontal locations of the 
domain. It is clearly seen from this figure that, indeed, the 
atmosphere keeps oscillating during all the simulation time, 
but the stationary regime is not reached since the variations 
in time are not strictly periodic, as was for the case of har¬ 
monic driving. There are few strong initial wavefronts that 
can be appreciated in the corona for t < 1000 s. Then, 
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the amplitude of the coronal oscillations becomes gradu¬ 
ally smaller. The dominant waves propagating in the corona 
are acoustic-like slow waves, and those with largest ampli¬ 
tudes are observed inside the vertical flux tubes (upper left 
panel for x = 150 km), consistent with the simulations of 
periodic driving. No appreciable fast magnetic waves are 
present after about t = 1000 of the simulation time, due 
to the reasons explained above. Fast acoustic waves prop¬ 
agating upwards from the lower turning point can also be 
appreciated at the left panels of Fig. at the bottom part 
of the domain. 

Since the pulse simulation does not reach the stationary 
regime, we are not able to calculate and analyze the average 
wave fluxes, as was done for the periodic driver simulations. 


3.4.2. Frequency distribution 


gradient in the Alfven velocity generated by the strong 

density jump. _ 

Similar results were found by Rosenthal et al. (2002), 
whose work is one of the first studies of wave propagation 
by means of MHD simulations. These authors find that 
the fast magnetic waves are refracted downwards by 
inclined magnetic fields while those waves propagating 
along the field lines in almost vertical magnetic fields 
continue propagating upwards unaffected by the field. 
In a three-dimensional atmosphere, Alfven waves may 
start to domina te the magnetic energy t ransport into the 
corona (see e.g. Ivan Ballegooijen et al., 2011 Khomenko 


fc Cally[ 2Q12|) and the magnetic energy may become of 


greater importance. Non-linear effects may also influence 
the amount of energy transported by the waves into the 
corona. 


The pressure pulse generates waves in a wide frequency 
range. The length of the simulation time series allowed us to 
perform the Fourier analysis of the simulations in order to 
determine the dominant period of the oscillations. For that, 
we Fourier-transformed the vertical velocity oscillations at 
each (x, z) location of the simulation domain and calculated 
the frequency corresponding to the maximum power in the 
spectrum. These frequencies are shown in Figure [T^ 

Figure reveals that the dominant frequencies of os¬ 
cillations below the photosphere lie in the 3-4 mHz range, 
almost independently of the magnetic structure. In partic¬ 
ular, 3-4 mHz oscillations are the dominant ones at the 
photosphere, 2 : = 0 km. The distribution of the dominant 
frequencies in the chromosphere is highly dependent on the 
magnetic structure. Oscillations reaching the chromosphere 
and corona are mostly in the 5-7 mHz frequency regime, 
or higher. One observes that some 3-4 mHz power reaches 
the chromosphere at the locations of the inclined arcades. 
Lower frequency oscillations also reach the corona along the 
external part of the flux tubes where the magnetic field is 
more inclined. At the internal part of the flux tubes, oscilla¬ 
tions of 6 mHz dominate. Curiously, the highest frequencies 
up to 10 mHz are observed above the transition region at 
the top parts of the inclined arcades and around the null 
point. 

4. Discussion and conclusions 

In this work, we have presented 2D numerical simulations 
for the study of the behaviour of MHD waves as they 
propagate from layers below the photosphere to the corona 
(up to 10 Mm). Our magnetic configuration (vertical flux 
tubes separated by an arcade-shaped structure) can be con¬ 
sidered as representative of a quiet sun region, in terms of 
the weak field strength and the varying orientation. In ad¬ 
dition, a null point exists where the magnetic field vanishes. 

After studying the wave behaviour at each atmospheric 
layer and different magnetic field inclinations, we conclude 
that, in a two-dimensional atmosphere, waves pass through 
the TR more preferably along vertical field lines. As the 
corona is a \ow-/3 plasma we have mainly acoustic energy 
propagating into it. Part of the magnetic energy stays 
concentrated in the lower corona just above the arcades. 
The lack of magnetic energy found higher up in the 
corona can be explained by the large refraction that the 
fast magnetic waves suffer in the TR, due to the strong 


We have also obtained that a large amount of energy 
propagates back downwards and upwards again in the 
atmosphere. On the one hand, there is the strong refraction 
of magnetic waves near the TR mentioned above. On the 
other hand, the strong gradient of the Alfven velocity in 
the vicinity of the null point also gives rise to a strong 
refraction of the waves around it, capable of sending waves 
downwards when propagating near it. These phenomena, 
together with the conversion of these waves into downward 
acoustic waves at the [3 = 1 layer and the turning point of 
these acoustic waves below the photosphere caused by the 
temperature gradient, lead to a continuously oscillating 
atmosphere, even if an instantaneous force is exerted. It is 
interesting to note that no periodical drivers are needed 
for this continuous oscillation. It must also be remarked 
that only the energy propagating inside the vertical flux 
tubes escape from the domain. This p henomenon was also 
suggested by Rosenthal et al. (2002), who stated that at 
some moment the downward propagating waves, due to the 
refraction in the TR, must eventually be reflected upward 
again, resulting in a resonant cavity. 


A frequency-dependent behaviour has also been obtained 
(see figure 13). Five-minute oscillations are channeled into 
the chromosphere and corona at the edges of the flux 
tubes where the magnetic field is more inclined, while 
three-minute oscillations penetrate into the higher layers 
in the more vertical magnetic fields. High-frequency waves 
can pass through the TR into the corona outside the 
vertical flux tubes and through the arcade. High-frequency 
transverse magnetic fast waves are the dominant ones in 
that region, even though their contribution to the magnetic 
energy is very smal l. This last result is similar to that of 


Fedun et al. (2011a), who also find an efficient transmission 
of high frequency waves into the corona inside a flux tube. 
Thus, it seems that waves with varying frequency may 
reach the corona under different conditions and magnetic 
field configurations. 
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LC. Santamaria et aL: Wave propagation from the sub-photosphere to the corona 



Fig. 7. Same as Fig. but for the simulation run with a 
horizontal periodic driving of 300 s. The black solid line is 
the layer where the cut-off period is 300 s. The movie is 
vailable online. 
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Fig. 8. Same as Fig. but for the simulation run with a 
horizontal periodic driving of 300 s. 
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Fig. 9. Same as Fig. but for the simulation run with a 
horizontal periodic driving of 300 s. The units of the bar 
are [erg/s/cm^]. 
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Vertical Mean Acoustic Flux 



Vertical Mean Magnetic Flux 



Fig. 10. Same as Fig. [^but for the simulation run with a 
horizontal periodic driving of 300 s. 



Fig. 11. Same as Fig. but for the simulations with an 
instantaneous pressure pulse. The movie of the wave prop¬ 
agation is available in the online version of the paper. 
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Fig. 12. Same as Fig. but for the simulations with an 
instantaneous pressure pulse. 



Fig. 13. Spatial distribution of the dominant frequencies of 
the vertically propagating waves for simulations with pres¬ 
sure pulse. Dashed-dotted black line marks /d = 1 contours, 
white solid line marks the transition region. 

































































